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We present tight-binding molecular dynamics simulations of the structural modifications that 
result from the "thermal spike" that occurs during the passage of a heavy fast ion through a thin 
diamond or amorphous carbon layer, and the subsequent regrowth upon cooling. The thermal spike 
and cooling down are simulated by locally heating and then quenching a small region of carbon, 
surrounded either by diamond or by a mostly sp^ bonded amorphous carbon network. For the case of 
the thermal spike in diamond we find that if the "temperature" (kinetic energy of the atoms) at the 
center of the thermal spike is high enough, an amorphous carbon region containing a large fraction 
of threefold coordinated C atoms (sp^ bonded) remains within the diamond network after cooling. 
The structure of this amorphous layer depends very strongly on the "temperature" of heating and 
on the dimensions of the thermal spike. Scaling is found between curves of the dependence of the 
percentage of sp^ bonded atoms in the region of the thermal spike on the heating "temperature" for 
different volumes. When the thermal spike occurs in an initially amorphous sample the structure of 
the damaged region after cooling exhibits the above dependencies as well as being a function of the 
structure of the original amorphous carbon layers. 



PACS numbers: 61.43.Bn, 61.72.Ww, 71.23.-k 

I. INTRODUCTION 

Carbon can bond in different hybridization forms, in- 
cluding diamond- like (sp^) and graphite-like (sp^), re- 
sulting in C based materials with extremely different 
physical and chemical properties. After equilibrium con- 
ditions are achieved, disrupted sp^ bonds may recon- 
struct as the more stable sp^ bonds. The passage of swift 
ions through matter causes severe bond breakage, which 
is generally referred to as a "thermal-spike" that accom- 
panies the stopping of the ion in matter. For the case 
of diamond, the local bond disruption of the sp^ bonds, 
and subsequent reconstruction upon equilibration, may 
result in the formation of an sp^ rich region along the 
ion-track. Indeed, the passage of MeV or GeV heavy ions 
through diamond-like (sp^ rich) layers has been shown 
by high resolution AFM and STM scans to result in lo- 
cal graphitization^ . This is evidenced by the formation 
of protrusions at the ion-impact points (swelling) and by 
the formation of electrically conductive channels. Both 
these factors indicate that the preferred sp^ bonding has 
occurred along the pathways of the ions, in the other- 
wise sp^ rich material^. The width of the transformed 
material, due to the impact of GeV Au ions, has been 
deduced from the data an was found to be of the order 
of 85 nm. Fig. 1 illustrates the material transformation 
(density and electrical conductivity) measured by atomic 
force and scanning tunneling microscopy. 

The damage inflicted on diamond by sub-MeV ion- 
implantation is known to result in electrical conductive 
regions along the ion track and around the stopping 
points of the ions'^. The nature of these conductive re- 
gions, which have been shown to give rise to hopping 
conductivity, is controversial. Prins postulates that they 
mainly contain a mixture of vacancies and interstitial C 
atoms which eventually collapse to from vacancy clus- 




FIG. 1: Surface morphology (left) and local electrical con- 
ductivity (right) of the transformed diamond-like amorphous 
carbon due to ion impact (taken from ref.^). 

ters ('vacloids')'*. On the other hand, Prawer and Kalish 
have concluded, based on measurements of the depen- 
dence of the conductivity of damaged diamond on ion- 
damage-implantation dose that the ion-affected diamond 
is highly conductive due to the formation of graphite-like 
regions along the ion-tracks^. 

Computer simulations of the damage in diamond 
caused by the energetic displacement of a few C atoms 
resulting in damage to a small volume in diamond^ (of 
1.4 nm^^) show that the size and nature of these regions 
depend on the initial recoil energy imparted to the C 
atoms and on the sample temperature. In some cases 
the formation of clusters of sp^ bonded C atoms forming 
a six-fold rings has been observed in these simulations. 

In the present study we simulate the passage of a swift 
ion through a diamond or a diamond-like tetrahedrally 
bonded amorphous C layer {ta — C) by locally "heating" 
a section of the sample up to very high temperatures 
(comparable to the temperatures expected to prevail in 
the thermal spike). We compute the sp^ and sp^ frac- 
tions in the affected volume and follow their evolution 
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function of cooling rate and time. We find that, 
upon cooling, a region which contains both sp^ and sp^ 
bonded C forms along the track. The nature of the mate- 
rial in this regrown region depends on the initial heating 
of the track and on its cooling rate. Under some con- 
ditions, complete regrowth back to perfect diamond is 
found, whereas in other cases, a sp"^ rich region remains 
along the ion track. Similar results are found when a 
thermal spike is simulated in amorphous carbon which 
is then allowed to cool down under a range of different 
conditions. The relevance of these computational results, 
obviously limited to small samples and to relatively short 
cooling times, to the actual ion passage through C based 
materials is discussed. 

We start below by briefly reviewing the physical pro- 
cesses involved in the formation of the thermal spike and 
the typical parameters of such a spike initiated in di- 
amond by the passage of a swift heavy ion; we then 
describe the computations performed, discuss their rel- 
evance to ion-implantation experiments; present and dis- 
cuss our results. 



II. BACKGROUND 

The "thermal spike" which accompanies the passage 
of an ion through matter can be viewed as a short-time 
local melting of the ion damaged region followed by a 
rapid cooling. Quenching from the liquid phase has been 
widely used to create amorphous materials (including 
amorphous carbon) in computer simulations, using em- 
pirical potentials^, ab initio^^^^ and tight-binding molec- 
ular dynamics^^"^^. It has been shown that such a proce- 
dure produces an amorphous carbon similar to material 
produced experimentally by the deposition of a carbon 
film either by evaporation in an electron beam, by car- 
bon arc^^~^^ or by ion beam deposition with ion energies 
between several ten and several hundred eV^''"^^. 

In contrast to this large body of computer simulations 
of amorphous carbon, there have been only a few simu- 
lations of the material resulting from damaging diamond 
by ion-implantation. Specific simulations of damage in 
material induced by ion implantation can be divided into 
two categories: (i) the commonly used TRIM code^°'^^ 
and (ii) the Molecular Dynamics (MD) technique^^ . The 
TRIM (TRansport of Ions in Matter) code is a Monte 
Carlo program that calculates the trajectories of the pri- 
mary ion and of the recoiling host atoms involved in the 
irradiation event. The required input parameters are 
the ion type and energy and the properties of the host 
material i.e. its composition, density, and the displace- 
ment energy of the constituent atoms. TRIM provides 
good estimates of the damage and of the implant dis- 
tributions and it is widely used in planning the dopant 
profile in ion-implantation experiments. However, TRIM 
cannot account for the structure of the material (amor- 
phous vs. crystalline) nor model the dynamic annealing 
of the implantation affected material. In particular, it 



cannot predict the nature of the defects in diamond be- 
cause of the variety of possible bonding configurations of 
carbon that may induce the diamond to graphite transi- 
tion. A modification of TRIM has recently been devel- 
oped by Sha'anan^^ which enables following the evolu- 
tion, in space and time, of the damage cascade caused by 
the passage of a single ion through matter. The results of 
the application of this code to simulate the passage of an 
IMeV Xc ion through diamond show that extremely high 
temperatures (of the order of many tens of thousand of 
degrees) are induced during very short times (less than 1 
ps) along the ion track. As we will discuss below, these 
resemble the time/temperature regime modeled in our 
computations of the damage inflicted on diamond and 
ta — C hy heavy ion irradiation. 

An MD study of radiation damage caused by the pas- 
sage of a single carbon atom through diamond was per- 
formed by Wu and Fahy^"^, using the Tersoff potential. 
In that study the computations were restricted to the in- 
vestigation of the damage threshold energy necessary to 
displace a single C atom from its lattice site. Smith^** 
simulated the results of the bombardment of diamond 
using a similar computational method but limiting the 
analysis to the ejection mechanisms of C atoms from a 
diamond surface (carbon self-sputtering) . Marks et. al.^^ 
have investigated the growth of a — C layers. By comput- 
ing the mean square displacement of the atoms during the 
cooling of the melt and the time required for it to relax, 
the life-time of the thermal spike in diamond was found 
to be less than 1 picosecond, for ion impacts below 400 
eV. In these latter study ab initio MD techniques were 
used. 

Saada et aP carried out MD simulations, using the 
Tersoff potential, of a damage region embedded inside a 
diamond matrix created by the energetic displacement of 
carbon atoms with a kinetic energy of up to 416 eV (8 
times the displacement energy of carbon atoms in dia- 
mond) in the bulk of a diamond crystal. Such regions 
are likely to appear in diamond at the final stages of the 
implantation-related damage cascade. The conclusion of 
that study was that partial amorphization of a volume 
with a radius of 1-2 nm occurs around the stopping point 
of the recoiling ion. With an increasing number of suc- 
cessively displaced atoms, the damage volume and the 
number of threefold coordinated atoms was found to in- 
crease. The structure of the damaged regions (i.e. the 
ratio of sp"^ and sp"^ bonded C atoms) was found in ref.^ 
to depend on the sample temperature. 

These computer studies strengthen our understanding 
of the processes occurring inside or at the end of an 'ion 
track' in diamond. However, the transitions from sp^ 
to sp'^ bonds within thermal spike regions have not yet 
been investigated in detail. The purpose of the present 
study is to investigate transformations of diamond and 
diamond-like amorphous carbon under the extreme con- 
ditions which prevail during and after the passage of a 
heavy fast ion through the materials. 



3 



III. TIGHT-BINDING SIMULATION DETAILS 
A. General description of the calculations 

We present a simulations of the thermal spike which 
occurs during the passage of a heavy fast ion through ei- 
ther diamond or amorphous carbon by locally heating 
a slab of the material to very high temperatures and 
allowing it to rapidly cool down. A large increase in 
volume accompanies the sp^- to sp^ bond conversion^^, 
but the volume of the damaged region in the center of 
the diamond is restricted by the surrounding diamond 
(or amorphous carbon) lattice, therefore simulations are 
carried out under constant volume conditions. In order 
to set the temperature, a randomly directed velocity is 
assigned to each atom. The value of the velocity is ran- 
domly selected in accordance with a Maxwell distribu- 
tion for the desired temperature. The velocities can be 
rescaled to be related to the ambient temperatures. Pe- 
riodic boundary conditions are applied to the samples in 
all three directions. In order to describe the interactions 
between carbon atoms we use a tight- binding modeP^, 
which has been shown to be transferable to successfully 
simulate different carbon polytypes: diamond, graphite, 
fullerenes^* as well as disordered carbon structures such 
as liquid and amorphous carbon phases^^'^^. The MD 
step was 0.5 x 10"^^ s. 

The method used here to distinguish between sp'^ and 
sp'^ bonded atoms is based on determination of the co- 
ordination number of each atom. The radial distribution 
function (RDF) of all carbon samples exhibits a clear gap, 
centered at about 0.19 nm, separating the first and the 
second peak corresponding to the first and second neigh- 
bor shells. All atoms within the sphere of radius 0.19 nm 
are thus assumed to comprise the nearest neighborhood 
of a given atom. Therefore, the number of neighbors 
of each atom within a distance of 0.19 nm determines 
the coordination number. Atoms that have four near- 
est neighbors within this sphere are assumed to be sp'^ 
bonded, while atoms which have only three are assumed 
to be sp"^ bonded. 

The Atomic Visualization package (AViz)^^ was used 
extensively in all stages of this project. A visualization of 
our carbon samples with color coding for different atomic 
bonding helped indentify clusters of either sp"^ or sp^ co- 
ordinated atoms, rings, graphite-like planes. Color pic- 
tures of the results of the present study can be found on 
the web site^°. 




FIG. 2; Initial configuration of the samples; grey (purple) 
balls represent the frozen atoms, white (yellow) balls represent 
atoms that are free to move. (Purple and yellow axe shown 
online) 



atoms), 2 X 2x 8 (256 atoms), 2 x 3x 6 (288 atoms) 
and 3 x 3x 6 (432 atoms) diamond unit cells. 

The "thermal spike" was created by heating the central 
layers up to "temperatures" of 1.2-2.6 eV (14000-30000 
K). The 32 upper and 32 lower atoms of each sample 
(the upper and lower layers of a height of 3.55 A) were 
frozen, i.e. the motion of these atoms was forbidden. 
The initial geometry of the samples is shown in Fig. 2. 
The outermost frozen diamond layers had the effect on 
their neighboring hot atoms of forcing them to return 
to their initial positions i.e. to regrow epitaxially back 
to diamond. In this way a temperature gradient from 
the edges to the center of the sample was created. Once 
the liquid phase reached equilibrium (this process was 
controlled by monitoring the total energy of the system), 
the central layers were cooled to room temperature at 
a cooling rate of 10 K/fs. The height of the hot layers 
varied between 4 and 6 diamond unit cells and the width 
of the hot layers varied from 2 to 3 diamond unit cells. 
The full time of cooling was 1.5 ps (3000 MD steps). 
However, the central layers ceased to move after the 750 
first MD steps, so the lifetime of our "thermal spike" was 
0.3 ps. This value is similar to the estimates of Marks^^ 
and Saada^, 0.2ps and 0.21ps, respectively. 



B. Details of the sample preparation 

1. Damaged diamond sandwiched between two layers of 
diamond 

Six different samples with a density of 3.5 g/cc, initially 
arranged as a perfect diamond crystal, were prepared. 
Their sizes were 2x2x6 (192 atoms), 2x2x7 (224 



2. Damaged amorphous C region sandwiched between two 
layers of frozen ta — C 

Four initial samples of diamond-like amorphous carbon 
at density 3.5 g/cc were prepared by heating diamond up 
to 5000 K followed by cooling of the liquid phase to 300 
K during 5 ps with an exponential cooling profile. Each 
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ta ~ C amorphous sample contained 192 atoms with a 
geometry 7.1 Ax 7.1 Ax 21.3 A. The samples contained 
66±4 % of sp^ coordinated atoms. In order to achieve 
a larger percentage of sp^ coordinated atoms (above 70 
%), two additional samples were prepared using the the 
following procedure: the amorphous structure was pre- 
pared with a density of 4.0 g/cc by quenching a liquid 
phase down from 10000 K to 700 K with a cooling rate 
of about 500 K/ps, then the density was reduced from 4 
to 3.5 g/cc by uniform expansion, after which the tem- 
perature of the system was raised to 2500 K, and there- 
after the structure was rccoolcd to room tcmperaturc"'^^. 
The two structures obtained at a final density of 3.5 g/cc 
consisted of 84 % and 77.5 % of sp^ bonded atoms re- 
spectively. 

To simulate a "thermal spike" in the diamond-like 
amorphous carbon samples, the central region of each 
sample was heated to "temperatures" of 1.2-2.60 eV per 
atom (14000-30000 K). The upper and lower layers of 
a height of 3.55 A each, were kept frozen. For the six 
amorphous carbon samples the height of the hot layers 
was 14.2 A (128 atoms in average). A representation of 
a typical initial sample geometry is given by Figure 2. 
Dark atom at top and bottom indicate the frozen layers, 
the light central ones are free to move and when heated, 
there is a temperature gradient from the edges to the 
center of the samples. Following the heating stage the 
structures were cooled to room temperature at a cool- 
ing rate of 10 K/fs. The time of calculations and the 
lifetime of the "thermal spike" were the same those de- 
scribed above for the case of the damaged region between 
two layers of perfect diamond. 



C. Sample equilibration 

We checked carefully, by monitoring of the total en- 
ergy of the system, that the heating time was sufficiently 
long to enable the liquid structure to reach equilibrium 
(early attempts to prepare samples of amorphous carbon 
by quenching from a liquid carbon that was not in equi- 
librium led to inconclusive results). 



D. Relevance to experiment 

Our aim is to determine the nature of transitions be- 
tween different carbon bonds as the result of the passage 
of fast ions in samples of different initial structures. To 
do this, we need to examine the region of passage and its 
intermediate neighbourhood on an atomistic scale. How- 
ever, this microscopic region is contained in a larger car- 
bon matrix, and while the microscopic region corresponds 
directly to that in the experiment, the surrounding ma- 
trix may differ in scale. We address the issues of system 
size and temperature in relation to the parameters of ex- 
perimental samples. 



In the simulations, on one hand, a maximum of 432 
atoms can be contained in a single sample due to com- 
putational limitations. On the other hand, the periodic 
boundary conditions mean that the sample is in a sense 
infinite, and, in particular, has no surfaces. Thus the 
simulation samples are both smaller and larger than the 
laboratory experiment. We are investigating the local re- 
gion around the ion path (which has a linear extent of 7.2 
A along the track directon prior to repetition). The ge- 
ometry was chosen so that this region is well surrounded, 
in the direction perpendicular to the track, by diamond 
or amorphous carbon matrix, which is frozen in place. In 
particular, there are 7.4 A between replicas of the tracks, 
a region that is sufficient to avoid interaction since the 
range of the potential is only 3 A. In order to ensure that 
the samples are of sufficient size to enable reliable inter- 
pretations, a series of different sizes were used for each 
situation and the effect of the scaling of the system size 
was carefully investigated. Although this is not exactly 
equivalent to the finite size scaling''^ that is commonly 
applied to phase transitions and critical phenomena, and 
is closer to the scaling used by Rosenblum et al^^ in a 
study of thermal stress at diamond interfaces. The excel- 
lent agreement of both of these scaling approaches with 
results from laboratory experiments gives us confidence 
in the concept. The collapse, as a function of system size, 
and the smoothness of the scaling that will be demon- 
strated below confirms that the results can be reliably 
extrapolated even further if needed. 

Temperature is introduced into the simulations as ki- 
netic energy. In the laboratory, the region around the ion 
track is extremely hot for a very short time. The high 
heating "temperatures" of the simulation are thus large 
amounts of kinetic energy imported the heated atoms. 
Here, too, we have studied the scaling of "temperature" 
with system size. Not surprisingly, the actual "tempera- 
ture" needed to attain specific structures is a function of 
system size. Again, collapse of graphs to a single func- 
tional form will be demonstrated and provides validation 
of our approach. 



IV. RESULTS OF TIGHT-BINDING 
SIMULATIONS 

A. Damaged diamond sandwiched between two 
layers of intact diamond 

1. Influence of the heating "temperature" 

When the diamond lattice in the central hot layers was 
melted, sp^ and sp bonded atoms appeared in the liquid 
layers. Upon cooling some (or all) of these sp^ and sp 
bonded atoms transformed back to sp'^ coordinated C 
atoms. The final fraction of threefold and twofold coor- 
dinated atoms depends on the kinetic energy per atom 
("temperature") in the heating process. The higher the 
heating "temperature" of the atoms the larger was the ra- 
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T (eV/atom) 


Fourfold(%) 


Threefold (%) 


Twofold (%) 


1.55 eV 


100 








1.72 oV 


95±5 


5±5 





1.98 eV 


79±5 


19±5 


2±2 


2.41 eV 


74±5 


23±5 


3±2 


2.58 eV 


62±5 


34±5 


4±3 



TABLE I: Fraction of the four-, three-, and twofold coordi- 
nated atoms in the heated layers after cooling in the sample 
with 192 atoms. 
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FIG. 3: Recovery rate showing the evolution of the fraction of 

two-, three- and fourfold coordinated atoms in 4 layers (128 
atoms), following heating to 1.725 eV and subsequent cooling. 



tio of sp^/sp^, that prevailed after cooling. Samples ini- 
tially heated to a relatively low "temperature" returned 
to diamond, whereas for hotter samples some non dia- 
mond bonds remained after cooling. 

For example, let us consider the 2 x 2x 6 sample 
containing 192 atoms, 128 of which were free to move 
(i.e. were heated and subsequently cooled). The per- 
centages of four-, three-, and twofold coordinated atoms 
in the hot layers after cooling from different "tempera- 
tures" are presented in Table I. Whenever the heating 
"temperature" was lower than 1.72 eV per atom all sp^ 
and sp bonds disappeared in the process of cooling, and 
the final, cooled, structure returned to a perfect diamond 
lattice, as shown in Fig. 3 where the percentages of dif- 
ferent hybridizations are graphed as a function of cooling 
"temperature" . 

At higher initial heating "temperatures" the heated 
layers remained, after cooling, partially or entirely amor- 
phous. The fraction of three- and two-fold coordinated 
atoms in the structure are shown in Fig. 4 which shows 
the percentages of sp"^ bonded atoms during the cooling 
process from different initial heating "temperatures" . 

The volume of the region remaining amorphous after 
cooling was found to increase with the heating "temper- 
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FIG. 4: The fraction of threefold coordinated atoms in the hot 
layers, shown during the process of cooling from the initial 
"temperature" of 1.725 eV and 1.98 eV/atom. The sample 
contained 4 hot layers (128 atoms). 



ature". In Fig. 5 visualization of some samples cooled 
from different heating "temperatures" are presented. We 
see that in the sample where the heating "temperature" 
was high (1.98 eV/atom) almost an entire layer of atoms 
which were free to move has transformed to amorphous 
carbon with 19 % of sp'^ bonds, while in the sample ini- 
tially heated to 1.38 eV/atom only a small central region 
is amorphous and the layers which are adjacent to the 
frozen part of the sample return to their crystalline di- 
amond structure. The reason for the dependence of the 
volume that remains amorphous on heating "tempera- 
ture" is that at lower "temperatures" the frozen outer 
layers of diamond extend their influence to larger dis- 
tances, therefore a larger volume of the damaged region 
reconstructs to the diamond structure after cooling. A 
tendency to segregation of the threefold and the fourfold 
atoms was observed in the samples heated to high tem- 
peratures when the volume of the region that remained 
amorphous was sufficiently large (almost all non-frozen 
atoms). 

The radial distribution function (RDF) of the C atoms 
in the hot layers is presented in Fig. 6. The second 
peak is visible in the figure due to the layers which have 
reconstructed to diamond structure after cooling. How- 
ever, as was explained above, the higher the "heating 
temperature" the larger the volume of the region remain- 
ing amorphous. Therefore all RDF peaks related to the 
higher heating "temperature" (dashed line) are lower and 
broader. For high heating "temperatures" the overlap- 
ping peaks are a superposition of two peaks that corre- 
spond to sp'^ bonded and sp^ bonded atoms. 

The angular distribution function of the hot layers 
demonstrates similar behavior, see Fig. 7. The peaks 
are lower and broader when the heating "temperature" 
is higher, this corresponds to a higher fraction of three- 
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FIG. 5: Amorphous carbon located between two layers of 
diamond. The samples contained 192 atoms (2 x 2x 6), tem- 
perature of heating was 1.2 eV/atom (left), 1.38 eV/atom 
(center) and 1.98 eV/atom (right). Grey (purple) balls rep- 
resent fourfold coordinated atoms, white (yellow) balls repre- 
sent threefold coordinated atoms, black balls (which appear 
only in the right most sample) represent twofold coordinated 
atoms. Boundaries between frozen layers and layers of atoms 
which are free to move are denoted by arrows. (Purple and 
yellow are shown online) 
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FIG. 6: Radial distribution function g(r) of central layers of 
the sample containing 4 hot layers (2 x 2x 6) cooled from 
two different initial "temperatures" . 



and twofold coordinated atoms. The angular distribution 
peak related to T =1.98 eV is located at 112° (closer to 
the ideal 120° for sp^ bonded C atoms), while the peak 
of the curve related to 1.72 eV is located at 109° (typical 
for sp^ bonded atoms). 



I20r 




60 80 100 120 140 160 



bond angle, 9 

FIG. 7: Angular distribution function of central layers of the 
sample containing 4 hot layers (2 x 2x 6) cooled from two 
different initial "temperatures". 



Number of layers 


Fourfold(%) 


Threefold (%) 


Twofold (%) 


4 


95±5 


5±5 





5 


78±5 


21±5 


1 


6 


74±5 


26±5 


0-2 



TABLE II: Fraction of the four-, three-, and twofold coordi- 
nated atoms in the samples with different height of the hot 
layer heated to 1.72 eV per atom. 



2. Effects of the size of the hot layers. 

The heating "temperature" is not the only factor that 

influenced the structure of the damaged region. If the 
height of the layer with atoms that were free to move 
was increased from 4 to 6 unit cells without changing the 
width of samples, the percentage of three- and twofold 
coordinated atoms in the damaged layers, when heated 
to the same temperature also increased (see Table II). As 
was explained above, the frozen layers assist the interme- 
diate layers in their vicinity to return to diamond lattice 
structure. The central atoms in the samples with a large 
number of the hot layers feel the influence of the dia- 
mond frozen layers to a lower extent, and therefore the 
volume of the region remaining amorphous after cool- 
ing was large, with a high percentage of sp^ coordinated 
atoms. In the small samples the central layers needed to 
be heated to very high "temperatures" to diminish the 
influence of the outer frozen layers. The highest "tem- 
perature" at which the sample reconstructed its diamond 
lattice after cooling also depended on the height of the 
hot layers. This temperature decreased as the height of 
the hot layers increased (see Table III). 

Graphs of the dependence of the percentage of sp^ 
bonded atoms on the heating "temperature" are plotted 
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Number of layers 


T* 


4 


1.63 ±0.04 eV 


5 


1.51 ±0.04 eV 


6 


1.38 ±0.04 eV 



TABLE III: The upper "temperature" of heating at which the 
samples reconstruct their diamond structure (T*). 
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FIG. 9: The same curves shifted to T* . The solid line is a 
result of square polynomial fitting. 



B. Simulation of molten and cooled amorphous 
carbon located between two frozen layers of 
amorphous cetrbon. 



FIG. 8: The dependence of the percentage of sp^ bonded 
atoms on the heating temperature for different heights of the 
hot layers. 



in Fig. 8 for the samples with different heights and same 
width of the hot layers. A striking feature of these curves 
is that they have the same shape for the different heights 
of the hot layer. All the curves can be fitted to a unique 
function f{T — T*), where T* is the upper temperature 
of heating, at which the sample can return to diamond 
after cooling. In Fig. 9 the collapse of the curves to a 
single scaling function is shown. 

The samples with different width (from 32 to 72 atoms 
in each layer) at the same height of the hot layers also 
demonstrate finite size effects. The percentage of three- 
and twofold coordinated atoms in the central layers is 
larger in the wider samples at the same heating "temper- 
ature" , the "temperature" at which the sample returns 
to diamond, T* , is higher for the thinner samples (see 
Table IV). 



Width 


T* 


1.72 eV 


2.07 eV 


32 atoms 


1.63 ±0.04 eV 


5 ±5 


20 ±5 


48 atoms 


1.48 ±0.04 eV 


19 ±5 


32 ±5 


72 atoms 


1.29 ±0.04 eV 


22 ±5 


41 ±5 



TABLE IV: The upper "temperature" of heating at which the 
sample reconstructed its diamond structure (T*) and percent- 
age of sp^ coordinated atoms in the samples cooled from 1.72 
and 2.07 eV per atom. 



Now let us consider structure of a central amorphous 
carbon layer sandwiched between two layers of frozen 
amorphous carbon. After cooling, the central dam- 
aged region remained amorphous, but the sp^/sp^ ra- 
tio changed. The dependence of this generated structure 
on the heating "temperature" was similar to that of the 
previous stage of the simulation (amorphous carbon be- 
tween two layers of diamond). Namely, if the central 
atoms were heated to higher "temperatures", the frac- 
tion of the sp"^ bonded atoms increased (correspondingly, 
the fraction of the sp^ bonded atoms decreased). Some 
of these structures were highly inhomogeneous, in par- 
ticular when the initial sample contained a large fraction 
of fourfold bonded atoms and was heated to a relatively 
high "temperature". In Fig. 10 the visualization of some 
of these amorphous carbon samples is presented. We can 
see that the frozen layers of the sample with 77.5 % of 
sp^ bonded atoms initially, contained after heating and 
subsecjuent cooling 75 % of sp^ bonded atoms, while the 
central layers after heating to 2.33 eV and subsequent 
quenching contained only 39 % of the fourfold atoms. 

The radial (RDF) and the angle (ADF) distribution 
functions of the central layers of this amorphous car- 
bon sample after cooling from the "temperatures" of 1.55 
eV and 2.33 eV are plotted in Figs. 11 and 12 respec- 
tively. The peaks corresponding to the higher tempera- 
tures of heating were lower, broader, and slightly shifted 
toward the shorter, graphitelike, bond length and larger, 
graphitelike, bond angle relative to those at the lower 
temperatures. This indicates a higher fraction of the sp^ 
bonded atoms. 

It is obvious that the percentage of differently coor- 
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FIG. 10: Samples of amorphous carbon located between two 
layers oita — C amorphous caxbon: the initial sample (77.5 % 

of sp^) (left), and the samples after cooling from 2.07 eV (in 
the center), and from 2.33 eV (right). Grey (purple) balls rep- 
resent fourfold coordinated atoms, white (yellow) balls repre- 
sent threefold coordinated atoms, black balls (of which there 
are only four in the right most sample) represent twofold co- 
ordinated atoms. Boundaries between frozen layers and layers 
of atoms which are free to move are denoted by arrows. 



dinated atoms in the frozen layers has influence on the 
results of the central region. The samples initially con- 
taining a larger fraction of the sp^ bonded atoms, con- 
tained a larger fraction of the sp^ bonded atoms after 
cooling (at the same temperature of heating). Fig. 13 
shows the final percentage of sp^ coordinated atoms as 
a function of heating "temperature" plotted for samples 
with different percentage of sp^ and sp^ bonded atoms 
before heating of the central layers. The data related to 
the smaller initial fraction of sp^ atoms have a lower final 
fraction of the sp^. 



C. The band gap 

During the calculations many samples of amorphous 
carbon with different structures were generated. The 
width of the band gap was "measured" for most of them 
by computing of density of states (DOS) as function of 
energy. For example, the graph of density of states for 
two samples of amorphous carbon, containing 86 % of sp^ 
bonded atoms and 53 % of sp^ bonded atoms is plotted 
in Fig. 14. The bandgap of the sample with 86 % of sp^ 
bonded atoms was estimated visually from the graph to 
be 2.3 eV. The sample with 53 % of sp'^ bonded atoms 
does not indicate the presence of any band gap, this sam- 
ple is expected to be electrically conducting. The width 
of the band gap depends not only on the fraction of the 
differently coordinated atoms, but also on their individ- 
ual configuration, i.e on the local clustering of sp^ bonded 
atoms. Fig. 15 shows a plot of the width of the bandgap 
as a function of the fraction of sp^ coordinated atoms in 
various ta — C samples. 

There is a noticeable scattering of the points between 
60 and 75 % of sp^ bonded atoms. We suspect that in 



70 
60 
50 
40 
30 
20 
10 
0. 





— T=1.55 eV/atonr 
-- - T=2.33 eV/atoiT 


1 

1 A. 




in 







0.5 1 1.5 2 2.5 

distance (Angstrom) 



FIG. 11: Radial distribution function g{r) of central layers 
of the sample containing 77.5 % of sp'^ coordinated atoms 
initially (79 % of sp'^ in the central layers) cooled from two 
different "temperatures": after cooling from 1.55 eV the cen- 
tral layers contained 78 % of sp^ coordinated atoms, while 
after cooling from 2.33 eV the central layers contained 39 % 
of sp^ coordinated atoms. 
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FIG. 12: Angular distribution function of central layers of the 
sample containing 77.5 % of sp^ coordinated atoms initially 
(79 % of sp^ in the central layers) cooled from two different 
"temperatures" . 



this interval the fraction of sp^ bonded atoms increases 
sufficiently to begin to create clusters. If the fraction of 
the sp'^ atoms is smaller than 25 %, than only separated 
sp"^ bonded atoms arc embedded in the sp^ amorphous 
network. If the fraction of the sp^ atoms is larger than 40 
%, the sp^ clusters begin to connect one with another. 
The band gap disappears when the fraction of the sp^ 
bonded atoms reaches 45 %. 
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FIG. 13: The final percentage of sp^ coordinated atoms as 
a function of floating "temperature" plottcxl for samples with 
different percentage of sp'^ and sp^ bonded atoms before heat- 
ing of the central layers. Circles in solid line-38 % of sp^ in 
initial amorphous sample (36% in the central layers), dia- 
monds in dashed line-36% of sp^ in initial amorphous sample 
(31% in the central layers), squares in dash-dot line - 30% of 
sp^ in initial amorphous sample (27.5% in the central layers), 
triangles in dashed line - 22.5% of sp^ in initial amorphous 
sample (21% in the central layers) stars in solid line - 16% of 
sp^ in initial amorphous sample (14% in the central layers). 
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FIG. 14: The density of states of the samples of amorphous 
carbon, containing 86 % of sp'^ bonded atoms and 53 % of sp^ 
bonded atoms. (£ = is chosen arbitrary.) 



V. DISCUSSION AND COMPARISON WITH 
EXPERIMENT. 

The experimental results^'^"'^^ described above show 
that the the high local energy deposition of the ions into 



FIG. 15: The width of the band gap as a function of the 

pcrccutage of sp'^ bonded atoms. The error bars represent 
an unaccuracy in visual determination of the band gap from 
DOS. 

diamond or amorphous carbon along their paths causes 
a rcarragement of the carbon atoms and can lead to a 
transformation of the insulating diamond or amorphous 
diamond-like sp^ form of carbon into a conducting layer 
rich in sp^ bonded carbon atoms. In the simulations we 
modeled the transformation of highly ordered diamond 
structures into amorphous structures with the presence of 
graphitic like sp^-bonds. In the samples which were ini- 
tially amorphous, rearragcmcnt of the amorphous struc- 
ture and an increase of the fraction of sp^ bonded atoms 
has occurred. 

The experiments'^'"^"' show that the resistivity, both in 
diamond and in amorphous carbon networks, decreases 
with ion dose. In diamond there exists a critical damage 
density Dc below which the diamond restores this pris- 
tine structure after annealing"'"'. At D = D2 a. sharp drop 
in resistivity is observed and at higher doses graphitic 
conduction is measured in the irradiated volume. Kalish 
et al^^ investigated the possibility that the sharp drop 
in resistivity could be explained in terms of a percola- 
tion transition between conducting islands produced by 
the impinging ion, but found that the Rys D dependence 
was not sharp enough to be consistent with such a simple 
percolation theory. By contrast, Hoffman et al^^ in their 
experiment on IkeV Ar irradiated diamond proved that 
the amorphization is a very sudden process and at D2 
the diamond collapses to an amorphous but still highly 
insulating form similar to amorphous diamond {ta — C). 
Only at doses exceeding those required for amorphiza- 
tion, does the amorphous sp^ structure transform to an 
amorphous sp^ structure which is electrically conduct- 
ing. In our simulation at low "heating " temperatures 
that corresponds to a slower or lighter ion impact, the 
diamond also reconstructs this structure after cooling. 
At higher heating "temperatures" the structure remains 
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amorphous after cooling and the fraction of sp^ bonded 
atoms in the damaged region increases with the heating 
"temperatures". The faster the ion the larger the dam- 
aged region. In our simulation increasing of the heated 
volume leaded to increasing of sp^ fraction in this ion 
track and to increasing of conductivity. 

Experimental results^^ show that in contrast to dia- 
mond , for ta — C where no crystalline template exists, 
there is nothing to inhibit the sp^ to sp^ bond conver- 
sion, hence there is no reason to expect a critical dose 
at which a sharp increase in the electrical conductivity 
should start. Moreover, the decrease of the electrical re- 
sistivity in ta — C occurs at much lower implantation 
dose than in diamond^^. The results of electron energy- 
loss spectroscopy allows to estimate the fraction of sp'^ 
bonds, /, in ion implanted ta — C. At low ion doses the / 
increases linearly with ion dose, at a dose 10^'' ions/cm'^ 
,/ becomes less sensitive and increases only slightly with 
dose. In our calculation ta — C becomes conductor at 
much lower heating "temperatures" than diamond, this 
means that in contrast to ta — C, diamond needs an ad- 
ditional energy to break its crystal lattice and lead to 
amorphization upon cooling. 

The conductivity in our amorphous carbon samples de- 
duced from the shrinkage of the band gap increases with 
percentage of sp^ bonded atoms in the damaged region. 
There are some slight discrepancies with experimental 
studies of electronic structure of damaged region. For 
example, McCuUoch et aP^ found a band gap of 2.5 eV 
at 85 % of sp^ bonded atoms, whereas we find that a 
value of 3.3 eV corresponds the fraction of 84 % of the 
sp^ atoms. 

VI. SUMMARY 

We have simulated a damaged region occuring in di- 
amond and in diamond-like amorphous carbon by local 
heating of a slab in the sample of thickness comparable 
to that created by the "thermal spike" due to the pas- 
sage of an energetic ion followed by quenching back to 
room temperature. The "temperature" of the heating 
varied from 1.2 to 2.59 eV/atom (14000-30000 K). For 
the case of diamond we found that at a heating "tem- 
perature" higher than a characteristic temperature T* a 
small region in the center of the sample remains amor- 
phous. The size and the structure of the amorphous re- 
gion depends on the heating "temperature" ; the size and 
the fraction of three- and two-coordinated atoms in the 



region increase with the heating "temperature" . When 
we increased the number of the free-to-move layers, the 
fraction of three- and two-coordinated atoms in the cen- 
tral layers increased. The curves of the dependence of 
the fraction of the sp^ bonded atoms in the central lay- 
ers on the temperature of heating have similar shape for 
different sizes of the hot layers. 

The same procedure was applied to simulate a dam- 
aged region in amorphous carbon samples. The central 
layers of each sample were melted at a very high tempera- 
tures (1.29-2.33 eV/atom) and then cooled by keeping the 
outer layers frozen. The percentage of three- and twofold 
coordinated atoms in the obtained damaged region was 
found to increase with heating "temperature" and with 
the initial fraction of the threefold coordinated atoms in 
the sample. Most of the structures of amorphous car- 
bon generated in different stages of our simulation were 
highly inhomogeneous and clustering of sp"^ bonded C 
atoms has been observed. 

The band gap measured for all the amorphous car- 
bon samples became narrower when the population of 
sp^ bonded atoms increased. It vanished when the frac- 
tion of sp'^ bonded atoms reached a value of 45 %. Be- 
yond this fraction the damage affected volume becomes 
conductive. 

The results obtained in this study support the exper- 
imental conclusions that in the damaged region formed 
after the passage of a heavy fast ion either through di- 
amond or amorphous carbon, amorphization and some- 
times conversion of sp^ bonds to sp^ bonds occurs. An- 
other conclusion which follows from the simulations is 
that after passing of a faster (or heavier) ion (that cor- 
responds to higher temperature of the "thermal spike" ) , 
the damaged region is larger and more sp^ bonds are 
transformed to sp"^. Although the present simulations do 
not exactly reproduce the conditions of the experimen- 
tal "thermal spike" , qualitatively our results are in good 
agreement with experiments. 
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